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We have measured the spin structure functions gi and g2 of '^He in a double-spin experiment by 
inclusively scattering polarized electrons at energies ranging from 0.862 to 5.058 GeV off a polarized 
■^Ife target at a 15.5° scattering angle. Excitation energies covered the resonance and the onset of 
the deep inelastic regions. We have determined for the first time the evolution of TiiQ^) = 
(;i(x,Q2)da;, T2{Q^) = g2[x,Q^)dx and d2(Q^) = x\2gi{x,Q^) + 'ig2{x,Q^)]dx for the 
neutron in the range 0.1 GeV^ < < 0.9 GeV^ with good precision. ri(Q^) displays a smooth 
variation from high to low . The Burkhardt-Cottingham sum rule holds within uncertainties and 
d2 is non-zero over the measured range. 

PACS numbers: 25.30.-c,1155.Hx,11.55.Fv,12.38.Qk 



During the past twenty five years, our understand- 
ing of Quantum Chromodynamics (QCD) has advanced 
through the study of the spin structure of the nuclcon. 
Measurements of the nucleon spin structure functions gi 
and (72 in deep inelastic lepton scattering (DIS) were used 



to unravel the spin structure of the nucleon in terms of its 
constituents, quarks and gluons, and test QCD. Among 
the important results are the finding that only a small 
fraction (about 20%) of the nucleon spin is accounted for 
by the spin of quarks jjy] , and the test of the Bjorken sum 
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rule 0, a fundamental sum rule of QCD, to better than 
10%. To make this latter test possible and determine the 
quark contribution to the total spin, it was essential to 
calculate the corrections necessary to evolve the Bjorken 
sum rule and the first moment of gi 0] to values acces- 
sible experimentally. To this end the sum rule, originally 
derived in the limit Q'^ ^ oo using current algebra, was 
generalized using the technique of operator product ex- 
pansion (OPE) in QCD IJ, la la la- In this connection, it 
was also realized that the Bjorken sum rule was a special 
case of a more general sum rule known as the extended 
Gerasimov-Drell-Hearn (GDH) sum rule Q which spans 
the full range of momentum transfer from = to 
0^ ^ oo. 

At small Q^, after subtracting the elastic contribution, 
f i(Q2) ^ ri(g2) - ri(Q2)ciastic^ linked to the anoma- 
lous magnetic moment of the nucleon k by 



tions: 



ri(Q')= / giix,Q')dx 



where Xq coincides with the nucleon pion threshold. The 
first term in the right hand side of ^ corresponds to 
the original GDH sum rule prediction The next 

term has been evaluated by Ji et al. using a heavy- 
baryon chiral perturbation theory (HBxPT) 1^ and 
by Bernard et al. using a covariant chiral perturbation 
theory (xPT) |i3,E3- 

At large (Q^ > ^ hrn] (Q^^ expressed in 
terms of a twist expansion |l3L Il4| : 



9Q2l«2+4d2+4/2 



where ag is the dominant, leading twist contribution. It 
is determined, apart from QCD radiative corrections ||6j, 
by the triplet gA and octet as axial charges and the net 
quark spin contribution to the total nucleon spin. These 
axial charges are extracted from measurements of the 
neutron and hyperons weak decay measurements [l5l |. 
Here a2 is a second moment of the gi structure func- 
tion and arises from the target mass correction [l^ . The 
quantities ^2 and /2 are the twist-3 and the twist-4 re- 
duced matrix elements. These matrix elements contain 
non-trivial quark-gluon interactions beyond the parton 
model. A first attempt at extracting /2 has been carried 
out by Ji and Melnitchouk using the world data but 
with poor statistics below — 1 GeV^. In QCD, d2 
and /2 can be expressed as linear combinations of the in- 
duced color electric and magnetic polarizabilities xe and 
Xb [l^ll^ when a nucleon is polarized. The above twist 
expansion may be valid down to = 0.5 GeV^ if higher 
order terms are small. 

We define ^2 as the second moment of a particular 
combination of the measured gi and 172 structure func- 



x^[2gi{x,Q^) + 3g2{x,Q^)]dx (3) 



3 / X- 
'0 



92{x,Q')-gY'^{x,Q^) 



dx 



where g^^ , known as the Wandzura-Wilczek ^21 term, 
depends only on gi 



gV'^ix,Q') = -g^ix,Q') + 



-dy. (4) 



The quantity d2 reduces to a twist-3 matrix element at 
large where an OPE expansion becomes valid. 

The advantages of measuring higher moments of the 
spin structure functions are twofold; 1) the kinematical 
region experimentally accessible gives most of the con- 
tribution to these moments, and 2) the matrix elements 
in the OPE of these moments can be calculated using 
lattice QCD [l^. 

Most of the previous measurements of Fi, r2 and d2 
were performed at well above 1 GeV^ where the 
higher-twist contributions are small compared to the pre- 
cision of the experiments. However, a good precision test 
of OPE requires precision data of P" starting from of 
about 0.5 GeV^ where multiparton interactions are im- 
portant. From = to perhaps = 0.2 GeV^, the 
moments predicted by the sum rules (e.g. spin polariz- 
abilities etc..) can be calculated in chiral perturbation 
theory (xPT) 0, 0| and can be tested against experi- 
ments. We do not expect OPE or xPT to be valid in the 
complete range of Q^; however, with time, lattice QCD 
may bridge the gap between these two limits. 

Finally, the 32 structure function itself is predicted to 
obey the Burkhardt-Cottingham (BC) sum rule 



r2(g')= / g2ix,Q^) dx = o 



(5) 



which was derived from the dispersion relation and the 
asymptotic behavior of the corresponding Compton am- 
plitude '2^1. This sum rule is also expected to be valid 
at all and does not follow from the OPE. It is a 
super-convergence relation based on Regge asymptotics 
as discussed in the review paper by Jaffe |2l|. Many 
scenarios which could invalidate this sum rule have been 
discussed in the literature jU |3. Ell' However, this 
sum rule was confirmed in perturbative QCD at order 
Oc with & g2(x,Q^) structure function for a quark target 



|27| . Surprisingly the first precision measurement of g2 at 
SLAC [2^ at Q2 = 5 GeV^ but within a limited range of 
X has revealed a violation of this sum rule on the proton 
at the level of three standard deviations. In contrast, the 
neutron sum rule is poorly measured but consistent with 
zero at the one standard deviation. 
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In this paper, we present measurements of the spin 
structure functions gi and (72 of ^He and the determina- 
tion of r5'(Q2)^ T^iQ"^) and d^{Q^) for the neutron below 
Q2 = I QeV2. The data on (Q^) and d^iQ^) provide 
a test of the latest results in xPT and permit a better 
extraction of f2- The data on T'l^iQ'^) allows us to make 
a more precise test of the neutron BC sum rule at < 1 
GeV2. 

The features of JLab experiment E94-010, from which 
these results are derived, were discussed in a previous 
determination of cttt' j the transverse-transverse virtual 
photoabsorption cross section 0|. We measured the 
inclusive scattering of longitudinally polarized electrons 
from a polarized '^He target in Hall A at the Thomas Jef- 
ferson National Accelerator Facility (JLab). Data were 
collected at six incident beam energies: 5.058, 4.239, 
3.382, 2.581, 1.718, and 0.862 GeV, all at a nominal 
scattering angle of 15.5°. The measurements covered val- 
ues of the invariant mass W from the quasielastic peak 
(not discussed in this paper), through the resonance re- 
gion and continuum. Data were taken for both longi- 
tudinal and transverse target polarization orientations. 
Both spin asymmetries and absolute cross sections were 
measured. More details can be found in 

In the Born approximation 51 (x, Q^) and (72(2;, Q^) are 
evaluated by combining data taken with opposite electron 
beam helicity and parallel or perpendicular target spin 
with respect to the beam direction. 



9i 

92 



MQ^iy E 1 
4a2 'E'E + E' 

MQ^j^^ 1 
4a2 2E'{E + E' 



-Act 



tan(6'/2)ACT_L 

E + E' cos ( 



(6) 



E' sin ( 



■Actj_ 



where Act||(^) = d^a^^^dndE' - d'^a'^^'^^^/d^dE' is 
the difference of cross sections for the case in which the 
target spin is aligned parallel (left-perpendicular) to the 
beam momentum. Here Ue is the electromagnetic cou- 
pling constant, 9 is the scattering angle, M is the nu- 
cleon mass, v is the transferred energy and E and E' are 
the initial and final energies of the incident and scattered 
electron respectively. 

The results of gi (circles) and (72 (squares) for ^He are 
shown in Fig. as a function of x for six values of 
in the range 0.1 GeV^ < < 0.9 GeV^. These struc- 
ture functions were evaluated at constant from those 
measured at fixed incident beam energies and angle by 
interpolation (filled symbols) , and for a few points (open 
symbols), extrapolation. The error bars represent the 
uncertainty due to statistics only, and the grey bands in- 
dicate the uncertainty due to systematic errors. The sys- 
tematic errors result from relative uncertainties of about 
5% in the absolute cross sections, 4% in the target polar- 
ization, 4% in the beam polarization, and 20% in the ra- 
diative corrections at every but the lowest incident beam 
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FIG. 1: gi (dark filled circles) and 32 (grey filled squares) of 
^He are plotted as a function of the Bjorken variable x for six 
values of . The points shown with filled (open) circles or 
squares were determined by interpolation (extrapolation). 



energy where it is 40%. They also include a contribution 
from interpolation and extrapolation. 

We notice a peak in gi and 32 due to the A1232 reso- 
nance, which decreases in magnitude with increasing . 
In the vicinity of that resonance gi and g2 are of almost 
equal amplitude and opposite sign. This is consistent 
with the fact that the A(1232) is an Ml resonance, lead- 
ing us to expect cflt' o: {gi + (72)13 to be highly su- 
pressed. We note also that in a region dominated by 
the coherent behavior of quarks and gluons (constituent 
quarks instead of current quarks) the Wandzura-Wilczek 
relation derived in DIS still holds, perhaps pointing to 
the role of quark-hadron duality in 171 . 

The integral Ti{Q^) for ■^He was computed for each 
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value of using limits of integration extending from 
the nucleon pion threshold to a value of x correspond- 
ing to = 2.0 GeV. To extract V^iQ'^) we followed the 
prescription suggested by Ciofi degli Atti and Scopetta 
in |3l|. where it is found, within the Impulse Approx- 
imation, that nuclear effects are quite significant when 
extracting the spin structure functions, but they reduce 
to at most 10 % in the extraction of r"(Q^). 
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FIG. 2: Results of (top panel), (middle panel), 
(bottom panels ) along with the world data from DIS and 
theoretical predictions (see text). 

The measured values of f"((5^) (solid squares) are 
shown in Fig. |21 (top panel) at six values of al ong 
with the world data (open symbols) from SLAC |35j 
and Hermes |36j . To obtain the total result of r^'(Q2) 
(solid circles ) an estimated strength in the DIS range 
4 GeV^ < < 1000 GeV^ using the parametriza- 
tion was added to the data obtained in the mea- 

sured region (solid squares). The size of each symbol 
indicates the statistical uncertainty while the system- 
atic uncertainties are illustrated with the dark grey band 



along the horizontal axis and include the uncertainty of 
the DIS contribution. We point out that at > 1 
GeV^ the elastic contribution is negligible and therefore 

f(g2) = r(g2). 

At low we show several xPT calculations; by 
Bernard et al. |ll[ without vector mesons (dotted line), 
by Bernard et al. [l3| including vector mesons and A de- 
grees of freedom (grey band) and by Ji et al. using the 
HByPT (dashed line). In the calculation of Bernard et 
al. [131 the grey band shows a range of results due to the 
uncertainty in the N7A form factor. The latter calcula- 
tion overlaps with a data point at = 0.1 GcV^. The 
GDH prediction is the slope depicted by the dot-dashed 
line at low Q^. 

At moderate and large we show the MAID cal- 
culation (3^1 (used to evaluate the resonance contribu- 
tion) combined with the DIS estimate from Bianchi and 
Thomas (solid line). The other calculations shown 
are by Soffer and Terayev [s^l (short-dashed line) and by 
Burkert and loffe [s^ (double-dashed line). 

While the MAID and Soffer and Terayev calculations 
are disfavoured, the result from the Burkert and loffe 
calculation agrees well with the data. More importantly, 
our data above — 0.5 GeV^ combined with the world 
data and future planned measurements [iol |4]| in the 
range 1 GeV^ < < 5 GeV^ will permit, in the frame- 
work of QCD and with better experimental constraints, 
to repeat the extraction of f2{Q^) higher- twist matrix 
element performed by Ji and Melnitchouk. 

We plot in Fig. [21 (middle panel) in the measured 
region (solid circles) and after adding the elastic con- 
tribution (open circles) evaluated using the Mergell et 
al. 01 parameterization of and G^. The open dia- 
monds correspond to the results obtained after adding to 
the open circles an estimated DIS contribution assu ming 



using the same method as described in [2 



92 = 92 

Nuclear corrections were performed following a proce- 
dure similar to that used in extracting F" including Q"^ 
dependent effects 0|- The solid line is the resonance 
contribution evaluated using MAID. The positive light 
grey band corresponds to the total experimental system- 
atic errors. The negative dark band is our best estimate 
of the systematic error for the low x extrapolation. 

Our neutron results (open diamonds) are consistent 
with the BC sum rule to within 1.7 standard deviations 
over the measured range. The SLAC E155x collabo- 
ration |0| previously reported a neutron result at high 
(open square), where the elastic contribution is neg- 
ligible, consistent with zero but with a rather large error 
bar. On the other hand, they reported a proton result 
which deviates from the BC sum rule by 2.8 standard 
deviations. Their quoted error bar in this case is 3 times 
smaller than that of the neutron result but still large. 

In Fig. El (bottom panel), d2{Q^) = d2{Q^) - 
^clastic j'Q2 -J shown at several values of Q^. The results 
of this experiment are the solid circles. The grey band 
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represents their corresponding systematic uncertainty. 
The SLAC E155x neutron result (open square) is also 
shown. The solid line is the MAID calculation containing 
only the resonance contribution. At low the HBxPT 
calculation |28| (dashed line) and the covariant xPT (dot- 
ted line) are shown. The two latter calculations overlap 
in the region shown. Furthermore, calculations of the 
covariant xPT with the vector mesons contribution (dot- 
dashed line) and the A degrees of freedom (long dashed 
line) are reported ^3 but are too close to the former 
XPT curves to be clearly seen at this scale. 

The lattice prediction ^ at Q"^ = 5 GeV^ for the 
neutron d2 matrix element is negative but close to zero. 
We note that all models (not shown at this scale) predict 
a negative or zero value at large Q^. At moderate 
our data show a positive ^2 indicate a slow decrease 
with Q^. 

In conclusion, we have made the first measurement of 
r'i(Q^), ^2iQ^) and d^(Q2) of the neutron from a 
regime where a twist expansion analysis is appropriate 
to a regime where xPT can be tested. The BC sum 
rule for the neutron is verified within errors in the inter- 
mediate range of due to a cancellation between the 
resonance and the elastic contributions. Our results 
suggest a positive contribution of the resonance region 
up to = 1 GeV^ of a size comparable to the SLAC 
E155x result. With time we expect our data, combined 
with future measurements to provide a challenging test 
of increasingly precise lattice QCD predictions. 
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